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Abstract 

We study possible contributions of heavy sterile neutrinos Vh to the 
decays r~ — > e^(/i^)7r^7r~. From the experimental upper bounds on their 
rates we derive new constraints on the — Vt mixing in the mass region 
140.5 MeV < rriy^ < 1637 MeV. We discuss cosmological and astrophysical 
status of I'h in this mass region and compare our constraints with those 
recently derived by the NOMAD collaboration. 
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1 Introduction 



Hypothetical sterile neutrinos Ug, blind to the electroweak interactions, were in- 
voked into particle physics in various contexts as possible means to resolve ob- 
served anomalies. The most prominent case deals with neutrino anomalies: solar 
neutrino deficit, atmospheric neutrino anomaly and results of the LSND neutrino 
experiment. Simultaneous explanation of all these three anomalies in terms of 
neutrino oscillations requires at least one sterile neutrino i^s In this case 

together with i^e,fi,T allows one to introduce three independent mass square 
differences associated with these three anomalies. Although the recent Super 
Kamiokande global analysis [0 disfavours active-sterile neutrino oscillations as 
a dominant channel this hypothesis is not ruled out. Along this line one may 
reasonably admit the existence of more than one sterile neutrino, say, one per 
generation, as it is suggested by various extensions of the SM. In this situation 
among the neutrino mass eigenstates, which are superpositions of the active z^e.^t.r 
and sterile z/s(j) weak eigenstates, one can encounter not only light neutrinos, but 
also heavy states u^. The question of viabihty of this scenario is the subject of 
experimental searches as well as cosmological and astrophysical constrains. 

Recently some phenomenological, cosmological and astrophysical issues of the 
intermediate mass neutrinos in the MeV mass region have been addressed 
in the literature 0]. This was stimulated by the attempts of explanation of the 
KARMEN anomaly @] in terms of a neutrino state with a mass of 33.9 MeV mixed 
with Ur . Although recent data of this collaboration ||^ have not confirmed this 
anomaly, the question of existence of heavy sterile neutrinos Uh remains open. 

These sterile neutrinos can be searched for as peaks in differential rates of 
various processes and by direct production of followed by their decays in the 
detector (for summary see Ref. [0 p. 361.). The can also give rise to signif- 
icant enhancement of the total rate of certain processes if their masses happen 
to be located in an appropriate region This effect would be especially pro- 
nounced in reactions that are forbidden in the SM. The lepton number/flavor 
violating (LNV/LFV) processes belong to this category. Many of them are strin- 
gently restricted by experiment and allow one to derive stringent limits on the Uh 
contribution. 

In the present paper we study the contribution to the LFV and LNV r- 
decays: r~ e~(yU~)7r+7r^ and r~ —>■ e^(yU"^)7r~7r~. The first (LFV) process 
can receive contribution both from Dirac and Majorana neutrinos while to the 
second (LNV) process only Majorana neutrinos can contribute. These processes 
are capable to provide us with unique information on the i>h — Vt mixing matrix 
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element Urh- In the mass region 140.5GeV < < 1637MeV, which we are 
going to study, the contribution to the considered r-decays gains resonant 
enhancement 0. This effect makes discussed r-decays very sensitive to presence 
of Vh with the masses m/j in the resonant region. Under certain assumptions we 
extract from the experimental data on these r-decays new constraints on the 
Urh matrix element in the resonant region for both Majorana and Dirac heavy 
sterile neutrinos v^- Up to recently there were no credible laboratory limits on this 
quantity in the MeV mass region [0 . Only recent data of NOMAD collaboration 
[ |T0[] allowed establishing constraints on Urh in the MeV mass region. In that part 



of the resonant region, which overlaps with the region probed by NOMAD, our 
constraints are more stringent by around two orders of magnitude. 



2 Neutrino mass matrix and interactions 

Consider an extension of the SM with the three left-handed weak doublet neutri- 
nos i/^j = z/^^, u'^^) and n species of the SM singlet right-handed neutrinos 
^'ri — i^Rii ■■■^Rn)- The general mass term for this set of fields can be written as 



1^ . „ 1 




-u'M'~''^u"' + ll.c. = --(u',iy'-){ t +H.C. 



1^ 



^'^m^iV^iVi + Yl.c. (1) 



=1 



Here M.l,M.r are 3x3 and n x n symmetric Majorana mass matrices, 

is 3 X n Dirac type matrix. Rotating the neutrino mass matrix by the unitary 

transformation to the diagonal form 

U^M^^'^U = Diag{m^i} (2) 

we end up with n + 3 Majorana neutrinos Ui = f/^jZ/^, with the masses mi,i. In 
special cases there may appear among them pairs with masses degenerate in 
absolute values. Each of these pairs can be collected into a Dirac neutrino field. 
This situation corresponds to conservation of certain lepton numbers assigned to 
these Dirac fields. 

The considered generic model must contain at least three observable light 
neutrinos while the other states may be of arbitrary mass. In particular, they 
may include hundred MeV neutrinos Uh, which we are going to consider in the next 
section. Presence or absence of these neutrino states is a question for experimental 
searches. 
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In this scenario neutrino mass eigenstates have in general non-diagonal LFV 
couplings to the Z-boson 

n+3 n+3 

KI^Pl^'^ = Z^ Ua^U:j^7]J,PLU, ^ J] V^J^l,PL'y^, (3) 

a=e,ii,T a=e,ii,T ij=l i,j=l 

where the last two expressions are written in the mass eigenstate basis. For the 
case of only three massive neutrinos one has Vmn = ^mn as a consequence of 
unitarity of In general Vmn is not a diagonal matrix and flavor changing 

neutral currents in the neutrino sector become possible at tree level. 

As is known the invisible Z-boson width is an efficient counter of the 
number of neutrinos. In the introduced scenario, despite possible presence of 
an arbitrary number of light neutrinos, the predicted value of the F^^ is always 
consistent with its experimental value. Indeed, if all the neutrinos are light with 
masses m^^j << Mz their contribution to the invisible Z-boson width can be 
written as 

n+3 

rL= E l%frf^ = rf^ E 5«A/. = 3rf^. (4) 

i,j=l a,f3=e,fj,,T 

With the SM prediction for the partial Z-decay width to a pair of light neutrinos 
YSM _ ig7.24±0.08 MeV this formula always reproduces the experimental value 
of Ff„^ = 498.8 ± 1.5 MeV. The chain of equalities in Eq. §) follows again from 
the unitarity of Thus, independently of the number of light neutrinos with 
masses rriu << Mz the factor 3 in the last step counts the number of weak doublet 
neutrinos. This conclusion is changed in the presence of heavy neutrinos with 
masses > Mz/2 which do not contribute to Fj„^. In this case the unitarity 
condition is no longer valid and the factor 3 is changed to a smaller value. 

Having these arguments in mind we introduce in the next section neutrino 
states Uh with masses in the hundred MeV region. These states can be composed 
of sterile and active neutrino flavors as described in the present section. 



3 T ^ e {fi )7T^7T decay rates 

Neutrinos contribute to the LNV r-decay /+7r^7r~ according to the lowest 

order diagrams shown in Fig.l. The t/-contribution to the LFV r-decay —>■ 
l~TT'^'K~ is determined by the tree-level diagram similar to the diagram in Fig. 
1(a). The loop-diagram analogous to the diagram in Fig. 1(b) is absent in the 
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(b) 



Figure 1: The lowest order diagrams contributing to r — ^ /+7r+7r decay. 



latter case. The diagram in Fig. 1(b) requires knowledge of the vr-meson wave 
function at long distances. As yet it is poorly known and would introduce into 
the calculations uncontrollable uncertainties. In the present paper we concentrate 
on the resonant mass domain determined in (|16D, where the diagram in Fig. 1(a) 
absolutely dominates over that in Fig. 1(b). For this reason we neglect this 
diagram latter on. The hadronic structure parameters necessary for calculation 
of the tree- level contribution Fig. 1(a) involves only one parameter of hadronic 
structure, the pion decay constant /tt, accurately known from experiment. For 
the T~ — > /^TT^TT" decay with / = e, /i we derive the following decay rate formula 



r(r l^Ti 7T ) = c ds 



E 



UlkUrkfTiuk 



k s-^Zk 



mt 



2^M Ids (y^wIz^) )dv{Y. 



—] (5) 

^ — 777, / rrj^ 



V{t — > / TT^TT ) = C ds 



E 



UlkUrk 



s — m 



vk 



s gH 



m 



2'' 



(6) 



The unitary mixing matrix IJij relates v[ = UijUh weak z/' and mass u neutrino 
eigenstates. The numerical constant is c = {G jp / 32) {n)~^ f^ml\Vud\'^ , where = 
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93MeV and = 1777MeV is the r-lepton mass. We introduced the functions 



G\z) = c^\z)z'' X [{z - xjf - xl{z + - 1)' - ^li^ + 1)], 

H\z^, Z2) = (1 + X] -Zi- Z2){x] - Z^Z2) + (7) 



(1 + X^i){zi + Z2- Xl) - 2{ziZ2 + X^) 



where Xi = mi/rrir, X{x,y,z) = x^ + + z^ — 2xy — 2yz — 2xz, 4>''{z) 
z, x'l)X^^'^{z, x^., xf). The integration hmits in Eqs. (|^), are 



m 



{xl - x^){xl - 1) + y{l + xf + 2x1) - Z/' ± <P\y) 



with y = s/ml. 

Assuming that there exist only hght neutrinos with masses m,yi << rrir we 
can rewrite Eqs. (j^), (j^) approximately in the form 



r(r- 

where 



I TV TV 



\{mu)lr\ 



-A 



(0 
+ ) 



I 71 71 



-A 



(0 



(9) 



(0 



ds f dv s 
— H ' 

V 



9 ' 9 . 
mt mt 



J s-^ mi 



(10) 



and 



{m^)ir = J2 UikUrk m^k, {ml)ir = ^ikUrk mlk- 

k k 



From atmospheric and solar neutrino oscillation data, combined with the tritium 
beta decay endpoint, one can derive upper bounds on masses of all the three 
neutrinos |]ll| m^^^^r < 3 eV. Then we have conservatively ||T^ 



|(m,),i| < 9 eV, \{ml)ri\<27 



'12) 
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With these upper bounds we obtain from Eqs. a rough estimate 

R.-i.-, - '^'"'T!''"?:'''' ; 3-0(1.2) X 10-". (14) 



This is to be compared with the present experimental bounds on these branching 
ratios 

K+V) ^ ^•9(3-4) X 10-6, 7^^-P^_^ < 2.2(8.2) x IQ-^, (15) 

The comparison clearly shows that the light neutrino contributions to the pro- 
cesses T~ /^vr'^Tr" are far from being detected in the near future. On the other 
hand experimental observation of t~ l^n^ir' decays at larger rates would 
indicate some new physics beyond the SM, or the presence of an extra neutrino 
state Uh with the mass in the hundred MeV domain. 

Assume that there exist neutrinos Uh with masses rrih in the interval 



i 140.5MeV <mh< Vs+ ^ 1637MeV. (16) 

These neutrinos would give resonant contributions to the processes l^n^7T~ 
since the first term in Eqs. (^ and the expression in Eq. (^) have non-integrable 
singularities at s = m^. Apparently, in the resonant region ( |T6| ) one has to take 
into account the finite width r,yh of the neutrino Uh. This is accomplished by the 
substitution rrih ruh — {i/2)Vyh in Eqs. (||), (|). As we will show in the next 
section the total decay width of the neutrino is small V^h « ^uh- Therefore 
in the resonant domain (|16D the neutrino propagator in Eqs. (^, (P) has a very 
sharp maximum at s = m|. The second term in Eq. (^), being finite in the limit 
T^h = 0, can be neglected in this case. This allows us to rewrite Eqs. (^), (|^) in 
the resonant mass region with good precision in a simple form 



r-(^- ^ /±7r^^-) ~ c7rG'(^o) '""'^;"'''^'"'' , (17) 
with Zq = {nih/mrY. 



4 Heavy sterile neutrino decays 

Heavy sterile neutrinos Vh can decay in both charged (CC) and neutral (NC) 
current channels Vh hh^2'JM{qiq2) and ^ilh^iM{qq), where I = e, fi 
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and M are mesons specified below. These decays are induced by the Lagrangian 
terms 



C = ^Uu irPL^^ w; + -^i^ f/„,f/:, va'PLy, z,, (18) 

where I = e,/i,r and Ua/B is the neutrino mixing matrix defined in Eq. (|]). 
Despite the fact that the NC term is of second order in the mixing matrix the NC 
decay channels of the heavy sterile neutrinos Uh, as will be seen, are as important 
as the CC one. 

Calculating the total width Ti^h we divide the interval (|16[) into two parts 

y^^l40.5MeV< rrih < rrin' = 958MeV domain (I), (19) 
rrir,' = 958MeV < rrih < Vs+ ^ 1637MeV domain (II), 

where m^/ = 958 MeV is the mass of the isoscalar pseudo-scalar meson ?7'(958). 

In the domain (I) we obtain T^h directly, by calculating all the Uh decay 
channels to the leptons / and mesons M one by one, using the following meson 
matrix elements 

(7r-(p)|J7^75M|0) = tV2Up^, {K-{p)\s^^j,u\0) = ~tV2fKP^, (20) 
{p'{p)\dlf,l5u\0) = V2mpfpPf,e*^, {K*'{p)\s-fp'j5u\0) = mK*fK*P,,e*p, 
(7r°(p)|M7^75M|0) = -(7r°(p)|(J7^75rf|0) = if^p^, 

(^(p)|m7m75m|0) = {v{p)\d-ft^l5d\0) = --(r?(p)|s7^75s|0) = 

{p%p)\u-f,u\0) = -{p'{p)\d^AO)=mpfpe;, 
{uj{p)\u-^pu\0) = {uj{p)\d-fpd\0) = m^/^e*, 

where is a vector meson polarization 4-vector; meson decay constants and 
masses are U = 93 MeV, fx = 113 MeV, fp = 153 MeV, fx* = 224 MeV, 

= 138 MeV, = 139.6 MeV, = 494 MeV, = 547 MeV, = 770 
MeV, = 782 MeV, rriK* = 892 MeV. 

This channel-by-channel approach can not be extended beyond the threshold 
of r7'(985) meson production since its decay constant /^/ is unknown. Above 
the 77' (985) threshold, in the domain (II), we approximate the decay rate of 
Uh — > l{i')M by the decay rate of the inclusive process Vh K^)liQ2 without 
perturbative and nonperturbative QCD corrections. This leading-order approx- 
imation should work for m^h >> A ^ 200 MeV. At lower masses a more viable 
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approach would be to relate the semileptonic Uh decay rate by the dispersion 
relations to the imaginary parts of the W and Z self-energies n(s) in analogy to 
the approach applied in the literature for the r z/+ hadrons inclusive decay 
|l^. However for our rough estimations we do not need this more sophisticated 
treatment and will use the above mentioned leading-order approximation. 
The decay channels of the Dirac neutrino Uh in the domain (I) are 

e"e+i/e, e"/i+i/^, /i"e+i/e, fi'fi'^Uf,, 
(CC) : Uh — > { e-7i+, fi-7r+, K+ , fx~K+, (21) 

e-p+, /i-p+, e-K*+, fx'K*+. 

(NC) : z//, — ^ I n' (22) 

[ Ui 7], Ui p^, UiUJ, 

where z/j = 1/1,1/2,1^3 are the three conventional light neutrino mass eigenstates. 
Latter on we neglect their masses since m^i << nth. 

Partial decay widths for these channels are readily calculated using the hadronic 
matrix elements in Eq. (pOl) . The results can be summarized as 

T{u, IP) = \Ui,\'^fpm^Fp{m,yp) ^ |f//,f rf), 

r(z/, IV) = \Ui^,\'av^fWhMyi, w) = It/^.frf \ (23) 



T{v, vy') = lop^ fyo (1 - y'yo)\i + 2yU ^ |oprr°\ 



where 



0= E U^hU:, (24) 

In Eq. (11) we denoted P = {7r+,K+}, V = {p+,ir*+}, P^ = {tt^^v}, = 
{p'^,Lj}, yi = rrii/mh, a'" = 1, a"^ = 1/3, a = 2sin^6'iy — sin^ 9w + 1/4, b = 
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sin^ ^iy(sin^ ^VK ~ V^)- The kinematical functions are 

' dz 



Hp{x, y) = l2f -{z- y^){l + x'- z)\''\l, z, x^)\''\<d, y\ z), 
i7« = 12/^(.-2,f)(l-zrv'r^ 



Fp{x, y) = Ai/^(1, x^ y') [(1 + x'){l + x' - y') - 4x% 
Fy{x, y) = X'/\l, x^ y') [(1 - x^ + (1 + x')y') - 2y% 



where the integration hmits are zi = y'^, Z2 = {1 — x)^ and w = 4yf. 

In the mass domain (II) of Eq. (^) we approximate the CC and NC decay 
rates by the rates of the inclusive reactions 

(CC): z/, ^rM2, {NC): Uf,-^u,qq, (25) 

where / = {e,fi} and q = {u,d,s}. The corresponding leading order decay rate 
formulas are 

T{u, ^ l-q,q-,) = |f/,,p^m^(|K.|2 + IK,/) ^ |[/,,|2r('^), (26) 

r(z/, ^iy,qq) = |ry,|2-^_rri^[9 - 16sin2^iy(l - sin^^^)] = Ir/.^r^^^). 

Here we neglected small cd, dc contributions. 

The following comment on the total NC decay rate T^i^ is in order. Summa- 
tion over all the NC channels gives 

r^T = KHj:h\' = K, E u^kU;,j:u:,u,, = (27) 

i=l a,f3=e,ii,T i=l 

a,f3=e,ii,T a=e,fi,T 

Here we assumed that there are only four neutrino mass eigenstates 2^1,2,3, ^'/i 
and used unitarity of the mixing matrix Ua/3- At the last step we neglected the 
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subdominant term ~ U^f^. Thus, both CC and NC decay rates are proportional 
to ~ and, as we mentioned at the beginning of this section, a priori are 

equally important. 

Taking into account this fact we collect the partial decay rates into the total 
z//i-decay width 

^uh = \Ueh? (r^ri + '^'^uh) + W^^h? {j^^ail + r^,'^^) + \Urh?T''uh (28) 

with 

iv^r^ + rf ) + e{m,, - m,) E rf") + e{m, - m^or^'^^, (29) 

M 

(r(^) + ri^)) + eim,, - m,) ^ r^'^^'") + 9{m, - m.or^'^^) (30) 

MO 

where the summations run over M = vr"*", K^, p"*", K*~^ and M° = 7r°, 77, p^, uo. 

If neutrinos Uh are Majorana particles, i.e. Vh = then both l~X{AL = 
0) and i^h — »• l^X'^[AL = 2) decay channels are open. This results in multiplica- 
tion of the right hand side of Eq. (|29| ) by factor 2. 



uh 



r 



(^) 

vh 



5 Bounds on heavy sterile neutrino masses and 
mixings 



Substituting the total z//j-decay width from Eq. (p8D into the resonant formula 
(p!7[), we can derive, from the experimental bounds (|I5|), constraints on the 
neutrino mass ruh and the mixing matrix elements Uah- In general these con- 
straints represent a hardly readable 4-dimensional exclusion plot. However under 
certain simplifying assumptions one can infer more valuable information on the 
individual size of the mixing matrix elements. In this paper we are interested in 
the Urh matrix element which is not constrained in the literature in the Vh mass 
range ([I6D (see [0]). Recently the NOMAD collaboration [jlOl obtained constraints 



for ruh < 190 MeV which overlap with a small part of this range. Therefore, it 
would be interesting to infer individual constraints on this mixing matrix ele- 
ment, at least roughly. A reasonable simplifying assumption would be to take 
\UTh\ ~ \U^^j\ ~ \Uej\- Then from the experimental bounds (|15D we derive a 2- 
dimensional — \Urh\^ exclusion plot given in Fig. 2 for the case of Dirac Vh- 
Multiplying the Dirac case limiting curve in Fig. 2 by the factor 2 one obtains 
the exclusion plot for the case of Majorana v^. For comparison we also present in 
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Figure 2: Exclusion plots in the plane jf/r/iP — n^'h- Here Urh and rrih are the 
heavy neutrino Uh mixing matrix element to and its mass respectively. The 
shaded regions are excluded by the big bang nucleosynthesis, the duration of the 
SN 1987A neutrino burst 0, by the NOMAD collaboration ^ and by r ^ 



decay [present result]. 



Fig. 2 the bounds from big-bang nucleosynthesis, SN1987A |^ and the NOMAD 
bounds []T0[. The cosmological bounds are shortly discussed in the next section. 

Note that at certain points of the mass range ([16|) the limits in Fig. 2 are 
incompatible with our assumption that lUrhl ~ \U^j\ ~ \Uej\ since other experi- 
ments already established limits like |f/ehP; I^^m/iP < 10^® at ruh = 350,400 MeV. 
Thus our limits around these points have to be treated with caution. Let us 
note that all known constraints on \Uih\ have been obtained under certain 
simplifying assumptions such as those that we used above. 

However, there is a special case when limits on \Urh\ can be extracted without 
any simplifying assumptions. This case is realized when there are experimental 
upper bounds on the rates of three processes involving rr, r/i and re lepton pairs. 
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These could be, for instance, decays 

B'^ r^r^TT^, r^/i^TT^, r^e^vr^. (31) 

Intersection of their resonant regions is 1917 MeV < mh < 4620 MeV. In this 
domain the decay rates can be written schematically as 

■p _ |^t/i|^ -p _ lUrhl'^lUfj.hl'^ 

" a\U,h\^ + b\U^^\^ + c\UrH\'' " a|f/e,|2 + 6|f/^,|2 + c|f/.^r 





\Urh\ 


2 


\Ueh\ 


2 


a 


\Ueh\ 




\' + c\U^h\ 


2 



(32) 



From these Eqs. we find 

\Urh\^ =arre + bTrf, + Crrr- (33) 

Notice that all the decay rates appear with positive degree. Only in such a case 
we can set upper bound on \Urh\'^ having upper experimental bounds r,-; < Expi. 
This would not be the case, for instance, in the system Tr^, Tre, Tee- 

At present the described approach can not be realized on practice due to 
absence of necessary experimental data. However the required upper limits on 
the semileptonic Sc-decays can be, probably, obtained in future experiments at 
the B-factories. Then, derivation of limits on \Urh\ free of simplifying ad hoc 
assumptions will become possible. 



6 Heavy sterile neutrinos in astrophysics and 
cosmology 

It is well known that massive neutrinos may have important cosmological and 
astrophysical implications. They are expected to contribute to the mass density of 
the universe, participate in cosmic structure formation, big-bang nucleosynthesis, 
supernova explosions, imprint themselves in the cosmic microwave background 
etc. (for a review see, for instance, Ref. |T^). This implies certain constrains 
on the neutrino masses and mixings. Currently, for massive neutrinos in the 
mass region ([TBI) , the only available cosmological constraints arise from the mass 
density of the universe and cosmic structure formation. 

The contribution of stable massive neutrinos to the mass density of the uni- 
verse is described by the "Lee- Weinberg" Q^h"^ — rrii, curve. From the requirement 
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that the universe is not "overdosed" this leads to the two well know solutions 
rrii, < 40eV and rriu ^lOGeV which seem to exclude the domain Eq. (|16D . How- 
ever for unstable neutrinos the situation is different. They may decay early to 
light particles and, therefore, their total energy can be significantly "redshifted" 
down to the "overdosing" limit. Constraints on the neutrino life times Tj,, and 



masses rrij in this scenario are found in Ref. [|T5[. In the mass region (|T6D we 
have an order of magnitude estimate 

T^^ < (~ 10^'^)sec Mass Density limit (34) 

Decaying massive neutrinos may also have specific impact on the cosmic struc- 
ture formation introducing new stages in the evolution of the universe. After they 
decay into light relativistic particles the universe returns for a while from the 
matter to the radiation domination phase. This may change the resulting density 
fluctuation spectrum since the primordial fluctuations grow due to gravitation 
instability during the matter dominated stages. Comparison with observations 
leads to an upper bound on the neutrino life time |]T6|. In the mass region ( |I6| ) 
one finds roughly 

T,^. < (~ 10^)sec Structure Formation limit (35) 



On the other hand, from our formula (pHf ) we estimate 



few sec < T^., Theoretical limit. (36) 

Thus massive neutrinos with masses in the interval ( |T^ ) are not yet excluded by 
the known cosmological constraints (|3^), (|35|) . 

Big-bang nucleosynthesis(BBN) and the SN 1987A neutrino signal may pre- 
sumably lead to much more restrictive constraints |^ . The impact of heavy sterile 
neutrinos Uh on the BBN is twofold: via their contribution to the cosmic energy 
density and via the secondary light neutrinos z/e, z/^, produced in z//i-decays. 
The first issue results in faster expansion which enhances the neutron-to-proton 
ration n/p. The secondary light neutrinos can produce also an opposite effect. 
The increased number of z/g maintains neutron-proton thermal equilibrium longer 
diminishing the n/p-ratio at the freeze-out point. There are also effects related 
to the I'e energy spectrum distortion and energy injection to the electromagnetic 
radiation from the e+e~ pairs produced in f/j-decays 0. The observed duration 
of the neutrino signal from the SN 1987A stringently constraints the energy emit- 
ted from the core in the form of penetrating particles like neutrinos. The sterile 
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neutrinos could affect this energy-loss limits and, therefore, their mixings with 
the active neutrinos have to obey certain constraints. 

Analysis of the BBN and SN 1987A constraints on z//j parameters has been 
made in Ref. 0. We present the corresponding exclusion plots in Fig. 2 for 
comparison with our limits. Unfortunately, the analysis of Ref. is restricted 
to the mass range m/j < 200 MeV for the BBN and rrih < 100 MeV for the SN 
1987A case. From the general appearance of the excluded bends one may expect 
that the extension of the SN 1987A bend can overlap with the region excluded by 
the r-decay, producing very stringent combined constraint around |?7r/iP < 10"^. 
Even stronger constraint would appear in the domain where the SN 1987A bend 
simultaneously overlaps with the BBN excluded bend. In this case a combined 
constraint would be, probably, more stringent than |?7t-/jP < 10~^^. 



7 Conclusion 

We discussed a generic model with sterile neutrinos which can contain heavy Ui 
mass eigenstates. If among them there are Uh states with masses 140.5MeV < 
rrih < 1637MeV they can resonantly contribute to r~ — > e^(/i^)7r'^7r~ decays. 

We derived decay rates for these processes valid both outside and inside of 
the resonant rrih region (0). In the latter case knowledge of the total i//i-decay 
width r^,/! is required. We calculated the r^,/! in the whole resonant rrih region 
(|16D taking into account both charged and neutral current decay channels. 

The effect of resonant enhancement of the heavy sterile neutrino Uh contribu- 
tion allowed us to extract from the experimental upper bounds on the r-decay 
rates rather stringent constraints on z/^ — Uh mixing matrix element \Urh\ in the 
whole resonant region (plGl) . The corresponding exclusion plot is presented in 
Fig. 2 together with the recent constraints from NOMAD collaboration and 



cosmological constraints 0. The r-decay constraints were derived under cer- 
tain simplifying assumptions about Uh-mixing with the other flavors. These 
or similar assumptions are always required when zz/j-mixing with one selected 
flavor is extracted from those processes which have been studied in the literature. 
We propose the set of processes Be rryr, r/xTr, reyr, yU/ivr, eevr which would al- 
low one to avoid such ad hoc assumptions and extract individual constraints on 
lUrhl, \ U^h\, \ Ueh\- Finally, we discussed cosmological implications of the heavy 
sterile neutrinos Uh and pointed out that extension of the known big bang nu- 
cleosynthesis and the SN 1987A constraints to the resonant region (|T6|) could 
significantly strengthen our constraints on \Urh\ by combining them with these 
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cosmological constraints. 
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